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We propose and characterize solid-state photonic structures where light experiences an artificial 
gauge field. A non-trivial phase for photons tunneling between adjacent sites of a coupled cavity 
array can be obtained by inserting optically active materials in the structure or by inducing a 
suitable coupling of the propagation and polarization degrees of freedom. We also discuss the 
feasibility of observing strong gauge field effects in the optical spectra of realistic systems, including 
the Hofstadter butterfly spectrum. 
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I. INTRODUCTION 

The effect of an external magnetic field on the dynam- 
ics of charged particles underlies a number of intrigu- 
ing phenomena in very different contexts, ranging from 
magnetohydrodynamics in astro- and geo-physics to the 
fractional quantum Hall effect in solid state physics [1]. 
Upon quantization, the eigenstates for non-interacting 
particles in a uniform magnetic field form a simple equi- 
spaced ladder of highly degenerate Landau levels in free 
space, while the interplay with a periodic lattice poten- 
tial was predicted to give rise to fractal structures in the 
energy vs. magnetic flux plane, the so-called Hofstadter 
butterfly So far, experimental observation of such a 
fascinating structure in ordinary solids has been hindered 
by the extremely high value of the required magnetic field 
intensity [3 . 

In recent years, an intense theoretical activity has in- 
vestigated the possibility of generating artificial gauge 
fields for neutral atoms by taking advantage of the Berry 
phase [4] accumulated by an optically dressed atom which 
adiabatically performs a closed loop in real space [5l [6] : 
the nucleation of a few quantized vortices in a Bose- 
Einstein condensate under the effect of an artificial gauge 
field has been demonstrated in the pioneering experi- 
ment by Lin et al. [7 . The combination of a gauge 
field with atom-atom interactions is expected to give rise 
to strongly correlated atomic gases that closely remind 
quantum Hall liquids jS]. 

In the meanwhile, experimental advances in the gener- 
ation and manipulation of photon gases in semiconduc- 
tor devices have opened the way to the study of collective 
many-body effects in quantum fluids of light [9 . The first 
reports of Bose-Einstein condensation [10 have been re- 
cently followed by the demonstration of superfluid flow 
around defects [Tl] and the hydrodynamic nucleation of 
vortices and solitons [12]. 

In this work, we theoretically investigate photonic de- 
vices where the orbital motion of the photon experiences 
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an artificial gauge field. Previous work in this direction 
has considered arrays of coupled optical cavities confin- 
ing single atoms [13 , topological electromagnetic states 
in gyromagnetic photonic crystals [14] , and time-reversal 
symmetry breaking effects for microwaves in circuit- QED 
devices [15]. In contrast to these works, our scheme can 
be implemented with standard solid-state photonic tech- 
nology in the visible or infrared spectral range [T6^. In 
combination with the on-going research on strongly cor- 
related photon systems p/Tj, it is expected to open new 
perspectives in the study of non-equilibrium many-body 
physics under strong magnetic fields. 



The basic idea of our proposal consists of imposing a 
non-trivial tunneling phase to photons by taking advan- 
tage of the polarization degree of freedom: this phase can 
be generated by an optically active medium embedded 
in the structure or can have a geometric nature. Geo- 
metric phases have been demonstrated in a number of 
configurations for propagating light [18, 19 . Here, we 
extend the idea to the case where photons are confined 
in a two-dimensional lattice and the phase is acquired by 
an evanescent wave while tunneling between neighboring 
sites. Two classes of devices, which can be built using 
passive dielectric materials with a real refractive index, 
are specifically considered. 

The first configuration is illustrated in Sec. |TI| and is 
based on an array of optical cavities: this scheme is 
suitable for observing general gauge field effects on pho- 
tons in the non-interacting regime and shows the inter- 
esting possibility of scaling the structure to any wave- 
length region. The second configuration is presented in 
Sec. |III| and is based on a single planar microcavity with 
a periodic lateral patterning. This scheme appears as 
most promising in view of combining the artificial gauge 
field with strong optical nonlinear it ies, so as to enter the 
regime of strongly correlated photon gases. To complete 
the study, Sec. |IV| is devoted to the discussion of some ob- 
servable quantities that can be used to extract the physics 
of quantum particles in strong gauge fields from experi- 
mentally accessible optical spectra. 
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II. THE FIRST SCHEME: ARRAY OF DBR 
MICRO CAVITIES 

This first configuration is sketched in Fig. [TJa) and 
consists of a two-dimensional array of optical cavities. 
Even if a similar physics can be observed in a wide class 
of systems, for concreteness we shall focus our attention 
on the specific case of monolithic distributed Bragg reflec- 
tor (DBR) cavities [20]. Within each cavity, light propa- 
gates along four arms oriented along the x, y axes, each 
of which is terminated by a DBR mirror. A polariza- 
tion preserving, weakly reflecting mirror oriented at 45 
degrees with respect to the x^y axes is located at the 
center of each cavity and serves to mix light in the differ- 
ent arms: in this way, the photonic eigenmodes of each 
isolated cavity will be linear superpositions of two stand- 
ing waves along the x and y axes. In the following we 
shall focus our attention on a single optical mode per cav- 
ity, e.g. the one consisting of a symmetric superposition 
of two standing waves along the two axes with a given 
circular polarization, say cr+ [21 . In order to suppress 
unwanted mixing of the circularly polarized (j± states by 
the central, oblique mirror, suitable circularly birefrin- 
gent layers can be used to lift the degeneracy of cf± cavity 
modes [22 . A more detailed discussion on the structure 
of the photonic modes within each cavity is given in the 
Appendix. 

Coupling between neighboring cavities occurs via 
evanescent wave tunneling across the separating DBR 
mirrors. In order to generate the tunneling phase respon- 
sible for the artificial gauge field, two options can be en- 
visaged. The first one involves inserting a pair of linearly 
birefringent half-wave slabs within each DBR mirror sep- 
arating neighboring cavities. The optical axes of the two 
slabs have a relative rotation angle 9 around the propa- 
gation direction. Propagation through the first half-wave 
slab transforms the incident <j+ light into cf- one; prop- 
agation through the second slab brings polarization back 
to cr+ , yet with an additional Pancharatnam phase factor 
^2ie q£ geometric nature [19]. On the other hand, cr+ light 
tunneling through the mirror in the backward direction 
will acquire an opposite phase factor e~^*^. Unwanted 
cavity-like resonances at the linearly birefringent slabs in 
the spectral vicinity of the cavity mode of interest can be 
ruled out by a careful choice of the slab parameters. 

The second choice involves a single slab of optically 
active medium in place of the pair of linearly birefrin- 
gent slabs. The tunneling phase is in this case gen- 
erated by the optically active material, which imposes 
phases ±a; An(i/2c to circularly polarized photons trav- 
eling across it in opposite directions. Here, d is the thick- 
ness of the optically active slab and An is the difference 
between the refractive indices experienced by the two he- 
licity states. 

To verify the existence of the tunneling phase when 
either a pair of linearly birefringent layers or an opti- 
cally active layer is embedded within a DBR mirror, we 
performed transmission matrix calculations for the sim- 
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FIG. 1: (Color online) Scheme of the coupled DBR microcav- 
ity configuration to generate the artificial magnetic field, (a) 
Two-dimensional square lattice of cavities: in the sketch, four 
cavities containing an oblique mirror each are connected by 
phase elements embedded between DBR mirrors. Within the 
Landau gauge, the artificial magnetic field arises from the x- 
dependence of the tunneling phase in the y direction, (b) Sim- 
plest two-cavity set-up. (c) For the two-cavity set-up of panel 
(b), transmission spectrum (blue solid line) for cr+ incident 
light; relative phase cj) of the field in the two cavity layers (red 
dashed line). The outer (inner) DBRs contain 19 (10) periods 
of alternate layers of refractive index ui — 3.6 and n2 — 2.9 
and optical thickness ABr/4, where Asr = 2nc/ujBv^ oj^v being 
the Bragg frequency. The cavity layers are Asr/S thick and 
have ricav = 2.9. They also contain a circularly birefringent 
medium of thickness da/XBr = 0.21/7r with n^^ = 1.2, 2. The 
angle between the optical axes of the half-wave slabs is ar- 
bitrarily chosen to be 37r/22. The same effect is observed in 
the presence of an optically active slab embedded in between 
the inner DBRs. 



plest two-cavity configuration schematically illustrated in 
Fig. [ijb). As one can see in Fig. [ijc), transmission is 
maximum when the incident frequency is resonant with 
one of the two eigenmodes of the coupled cavity system. 
The presence of the non-trivial tunneling phase is ap- 
parent as a non-vanishing relative phase of the electric 
field in the two cavities: in standard configurations with- 
out any phase element, the two eigenmodes correspond 
to the symmetric and antisymmetric combinations of the 
isolated cavity modes, while now each of the eigenmodes 
exhibits an additional (j) relative phase between the two 
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cavities. 

In a tight-binding model, the tunnehng phase can be 
described in terms of a hopping Hamiltonian of the form 
He = —Je^^ ^\i^L +h.c., with J being a real and positive 
coefficient quantifying the strength of tunneling; Bl^r are 
the cavity mode operators for respectively the left and 
right cavities: it is straightforward to see that the eigen- 
modes of He are indeed linear superpositions of the bL^R 
isolated cavity modes with a <p phase difference between 
the two cavities. 

The form of this coupling term directly extends to the 
full two-dimensional geometry shown in Fig. [ija), where 
the tunneling phase between each pair of neighboring cav- 
ities can be independently tuned by the relative orienta- 
tion of the two half-wave slabs. As usual, the hopping 
phase between the neighboring z, j cavities (equal to (j) 
in the previous example) can be written in terms of an ar- 
tificial gauge potential A as = — | J^' A-dl, where e is 
the elementary charge and the integral is performed along 
the segment connecting the cavities. A no n- vanishing ar- 
tificial magnetic field then appears in the photon dynam- 
ics whenever the sum of the tunneling phases around a 
closed loop is non-zero (modulo 27r). 



III. THE SECOND SCHEME: LATERALLY 
PATTERNED PLANAR DBR MICROCAVITY 



The second scheme is based on a planar DBR micro- 
cavity architecture. Full three-dimensional confinement 
of the photon in micron-sized wells has been demon- 
strated by means of a lateral patterning of the cavity 
layer thickness [23] . A periodic array of photon boxes can 
be obtained by means of a two-dimensional periodic rep- 
etition of the elementary well. As usual in Hubbard-like 
models, a tighter confinement within each well allows to 
enhance photon-photon interactions: with a sub-micron 
confinement, a strongly interacting photon regime is ex- 
pected to be accessible using state-of-the-art semiconduc- 
tor technology [T71I24] . 

In addition, we assume a position-dependent vector 
field to be present that couples to the effective spin- 1/2 
system describing the photon polarization state in the 
a± basis (defined with respect to the cavity growth axis 
z): the z-component of this field can result from a static 
magnetic field that splits the a± polarization states by 
inducing a circular birefringence in the cavity material or 
by splitting the Zeeman components of an exciton state 
to which the photon is coupled [25 . On the other hand, 
the X and y components of the vector field (i.e. the ones 
that mix (j± polarization states) can be generated via a 
linear birefringence of the cavity material induced e.g. by 
a mechanical stress [26J or by a sub-wavelength grating 
imprinted on the cavity [27] . 

Specifically, we consider the one-dimensional configu- 
ration sketched in Fig. [2ja) that gives rise to a Hamilto- 
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FIG. 2: (Color online) Laterally patterned planar microcavity 
scheme to generate the artificial magnetic field, (a) Spatial 
dependence of the scalar potential Vsc{x) (solid line) and the 
amplitude of the xy component of the vector field (dashed 
lines) ; different dashed lines correspond to the components of 
amplitude Vs{x — Xj) centered at Xj having different angles 
with Cmax = 27r/3. Effective photon mass m was taken to 
be 5 X 10~^ times the electron mass. The z component of the 
vector field is Vz{x) — —1.21 meV. (b) Spatial dependence of 
the ground state expectation value of different spin compo- 
nents {Si) — \{cri). Inset shows the corresponding loop on 
the Poincare sphere. 



nian of the form 



H = 



2m 



(1) 

where m is the effective photon mass along the two- 
dimensional cavity plane, a's are the Pauli matrices in 
the two-dimensional spin space spanned by the a± po- 
larization states. = exp(— Z(j;^C/2) are the rotation 
operators around the z axis. 

The scalar potential Vsc{x) [solid line in Fig. [2|a)] 
stems from the lateral patterning of the cavity layer 
thickness and confines the photons in two square wells; 
the additional barrier in the center serves to cancel the 
localizing effect of the vector field in the xy plane. The 
component Vz{x) of the vector field coupling to the z 
component of the effective spin is assumed to be constant 
in space. 

The component of the vector field along the xy plane 
is localized in between the two wells and is modeled as a 
superposition of several {j = 1, 2, . . . , jmax) localized po- 
tentials of amplitude Vs{x — Xj) centered around neigh- 
boring positions Xj and oriented in different directions 
making angles Q with the x direction gradually varying 
from to a maximum value Cmax- 

The smooth variation of the angles Q is intended to 
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ensure adiabaticity: For large enough amplitudes of the 
vector field, the two spin states are energetically split. As 
a result, the photon polarization is able to adiabatically 
follow the local ground state determined by the direction 
of the local field and traces a closed loop on the Poincare 
sphere. On general Berry phase arguments [4 , we can 
anticipate that tunneling between the wells will involve 
a geometric phase which is equal to half the solid angle 
ft subtended by the closed loop. 

This expectation has been verified by a numerical cal- 
culation of the ground state of the Hamiltonian ([T]) by 
means of an imaginary-time evolution. The ground state 
is localized within the potential wells determined by the 
scalar potential and the local expectation value of the ef- 
fective spin operator indeed follows a closed loop on the 
Poincare sphere as shown in the inset of Fig. [2|b). As 
expected, the relative phase of the ground state wave- 
function in the two wells is found to be very close to the 
value (/) = Q/2 predicted by the adiabatic model. 

A full two-dimensional lattice of wells can be obtained 
by repeating this building block along both directions. 
With a suitable tuning of the hopping phase between 
pairs of neighboring wells, the photon turns out to expe- 
rience a non-trivial artificial gauge potential A. 

IV. OBSERVABLES 

After having discussed possible methods of creating 
an effective magnetic field for photons in a lattice, we 
will now turn to its observable consequences. For the 
sake of simplicity, we shall concentrate on the case of a 
uniform magnetic field and vanishing photon-photon in- 
teractions. In our calculations we consider a finite-size, 
two-dimensional square lattice within the tight-binding 
limit and we include the pumping and loss terms describ- 
ing the coupling of the cavity system with the outside 
world in terms of a master equation in the standard form 

dtp=i[p,H]/h^c[p] m- 

In the case of non-interacting photons, the tight- 
binding Hamiltonian of the isolated system has the fol- 
lowing single-particle form 

H = Y^ huJoblk-hJ blbj(^''^''^Yl [^^^(^) + ^•^•] ' 

(2) 

where h\ {hi) is the bosonic creation (annihilation) op- 
erator for site i. The hopping phase stems from the 
artificial gauge field; ujo is the natural cavity frequency 
and J is the tunneling strength between nearest neighbor 
sites. We assume the coherent driving term Fiit) to be 
monochromatic at frequency ujp and to act on the single 
site n, Fi(t) = F din e~*^p^ Photon losses at a rate 7 are 
included via the standard Lindblad term 

i 

As we are considering a non-interacting system, the state 



of the field is a product of coherent states on each site 
with an amplitude f3i = (hi) determined by the corre- 
sponding classical field evolution equations 

z A = (a;o - ^7/2) A - J ^ e*^- + F, (t) , (4) 

(j) 

where the sum over (j) is restricted to the nearest neigh- 
bors of site i. 

As usual in optical devices, the amplitude of the emit- 
ted light by each site is proportional to the bosonic oper- 
ator bi'. differently from the standard paradigm of quan- 
tum mechanics, the phase of the photonic wavefunction 
is then an experimentally accessible quantity, which is 
sensitive to the gauge potential A and not only to the 
magnetic field V x A. 
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FIG. 3: (Color online) (a) Average total number nr of trans- 
mitted photons as a function of The system 
considered is a 3 x 3 lattice with flux quanta per plaquette 
a = 1/3 and 7/J = 0.05. Different (green, blue, red) curves 
correspond to a driving of amplitude F/J = 0.005 localized 
on different (2,3,5) sites according to the enumeration given 
in the inset. Solid black line shows the sum over all nine 
sites, (b) Real-space profile \ip''^\r)\'^ of eigenmodes 1^^), 
where Gaussian real-space basis functions were assumed at 
each site for illustrative purposes. 

As a specific example, we consider the case of a uni- 
form magnetic field with the number of fiux quanta per 
plaquette a = (27r)~^ ^^j, where the sum is along a 
closed loop surrounding the plaquette. The total num- 
ber of photons riT = ^i(bih) present in the system is 
plotted in Fig. [sj^a) as a function of the pump frequency 
Aujp = UJp— UJo for a 3x3 lattice at a = 1/3 with hard- wall 
boundary conditions, as is relevant for experiments. In 
the laterally patterned microcavity configuration of Fig. 
|2j riT is proportional to the total transmitted intensity 
across the system. 

In Fig. IsFa), the different colored curves (^) cor- 
respond to^ocalized driving on different sites n. All 
curves exhibit peaks at frequencies corresponding to the 
eigenmodes of the system. For each eigenmode, the peak 
strength is proportional to the weight of the eigenmode 
on the driven site n. This physics is summarized by the 



5 



explicit expression 



^T,(n) 



^(6l6.)(„) = |F|2^ 



(0|2 



(5) 

where ipn"^ is the component on site n of the wavefunc- 
tion corresponding to the eigenmode / of the hopping 
Hamiltonian [i.e. the second term in Eq. and Ei is 
the corresponding eigenfrequency. Examples of the wave- 
functions V^^^^ of different eigenmodes are shown in Fig. 
ib). 

Analogously, the local field on site i under a localized 
drive on site n is given by 



^ Aojp + ij/2 - El ' 



(6) 



which can be recognized as F times the Green's function 
G(z, n, Aujp -\-ij/2) of the non-lossy system, measured at 
frequency Aujp + ^7/2 and positions z,n. Green's func- 
tion methods have been widely used in the literature to 
study, e.g., the effects of disorder on the Hofstadter spec- 
trum [29j . 



but finite losses. In Fig. [4]^ a), this quantity is plotted 
as a function of both AcUp and a for a 10 x 10 lattice. 
In addition to the clearly visible butterfly structure that 
closely resembles the infinite-size case, one can also rec- 
ognize a series of low intensity lines appearing within the 
largest energy gap. Direct inspection of the eigenstates 
shows that these lines correspond to edge states. From 
Fig. [4]^b), it is apparent that the separation between 
neighboring lines has an almost constant value Aa ap- 
proximately equal to I /A, where A= (L — 1)^ is the area 
enclosed by the outermost sites of an L x L lattice. This 
value of Aa corresponds to a change of magnetic flux 
across the whole lattice by one flux quantum and can be 
interpreted by the Streda formula for the quantized Hall 
conductance [30^. 

Another interesting feature of the Hofstadter physics is 
the spatially periodic structure of the ground state wave- 
function for rational values of a [31 . In our photonic 
system, this can be experimentally studied by tuning the 
pump frequency on resonance with the lowest frequency 
peak and collecting transmitted light from each site in a 
spatially-selective way. As an example, we show in Fig. 
|4jc) that for a = 1/3 the pattern exhibits a simple peri- 
odicity of 3 sites. 
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FIG. 4: (Color online) (a) Color plot of the total transmission 
riT as a function of flux quanta per plaquette a and pump 
frequency Aujp/J for a 10 x 10 lattice with F/J = 0.005 and 
7/J = 0.05. (b) cut of the color plot along the constant Aujp 
line indicated in (a) as a black line, (c) Photon occupation 
number (blbi) pattern for a pump resonant with the ground 
state and localized on a single site at the lattice center. 21x21 
lattice, magnetic flux quanta per plaquette a = 1/3, 7/J = 
0.01. 

To clearly illustrate the Hofstadter butterly, it is use- 
ful to plot the sum of transmission spectra over all pos- 
sible experimental realizations in which only one site is 
pumped at a time, as shown by the black line in Fig. 
[sja). It can be shown that in the limit of vanishing losses 
this quantity Xl^^T,(n) is proportional to the density of 
states of the Hofstadter spectrum (modified by finite-size 
effects), so it is natural to expect that the main features 
of the spectrum will be retained in the presence of weak 



V. CONCLUSIONS 

In this work, we have proposed two configurations 
where the photon experiences an artificial gauge poten- 
tial in a solid-state photonic device. While the first 
scheme appears suitable for a first experimental demon- 
stration of artificial gauge fields for non-interacting pho- 
tons, the second one is promising in view of combining 
the gauge field with strong photon-photon interactions. 
We have pointed out observable consequences of the arti- 
ficial gauge field in the experimentally accessible optical 
spectra of the device and have identified clear signatures 
of the Hofstadter physics at strong magnetic fields. Fu- 
ture theoretical work will aim at extending this study to 
the interacting regime where novel non-equilibrium fea- 
tures of strongly correlated quantum Hall fluids of light 
are expected to appear. 
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Appendix A: Structure of the photonic modes 
within each cavity 

In this Appendix, we provide more details on the struc- 
ture of the photonic eigenmodes within each cavity of the 
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two-dimensional array considered in Sec.jTTj In particular, 
we analyze the mechanism of the coupling between light 
beams propagating along the x and y directions. The 
element that couples the two directions is a partially re- 
flecting mirror inserted inside the cavity making an angle 
of 45 degrees with the x and y axes as shown in Fig. [5] 




FIG. 5: Schematic arrangement of mirrors forming the cavity. 
Central mirror is located at 45 degrees with respect to x and 
y axes, its center being at a distance L from the external 



We expand the electric field along each of the four arms 
inside the cavity as two counter-propagating plane- waves 
and use the appropriate scattering matrices to match the 
incoming and outgoing field amplitudes at the central 
and external mirrors. This matching yields the following 
matrix equations 



t r 
r t 



t r 
r t 
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where t (r) is the transmissivity (refiectivity) of the ex- 
ternal mirrors defined as the ratio of the transmitted (re- 
fiected) amplitude to the incident amplitude, to {vo) is 
the transmissivity (reflectivity) of the central mirror em- 
bedded in the cavity, and (p = uoLn/c is the phase that a 
wave gets after it travels a distance L in the cavity, n be- 
ing the refractive index of the cavity material. The field 
amplitudes are defined in the immediate vicinity of the 
external mirrors. Here, for simplicity, we assumed that 
the media inside and outside of the cavity are the same 
and further took all mirrors to be symmetric. This con- 
dition, together with the unitarity of scattering matrices 
(dictated by flux conservation) is satisfied by choosing r 
and To to be purely imaginary and t, to to be real. 

So, the mirrors can now be characterized by only one 
parameter, as transmissivities and refiectivities are re- 
lated through r = Wl — t'^ and To = iy^l — t1. The 
i factor in the reflectivities means that upon reflection 
from the mirrors the fields acquire a phase factor of e*^/^. 
The six matrix equations [Eqs. ( |A1|A6 )] written above 
provide us with twelve linear equations for the twelve un- 
known fields all scaled by the incident field Ei. In what 
follows we will adopt the convention that Ei = 1. 

Let us first suppose that the central mirror is perfectly 
transmitting, i.e. to = l,ro = 0. It is obvious in this 
case that no field will develop along the x direction (i.e. 
^3,4,7,8 = 0) as the incident field is propagating along y. 
Solving the equations one can find the transmitted field 
intensity to be \Et^\^ = |t|^/(l + \r\^ + 2|rp cos4(/p). In 
order for this quantity to be maximum so that we could 
say a cavity mode develops, the condition cosAcp = — 1 
should be satisfied yielding 



UJ = 



7r(2A/' + l)c 



(A7) 



where N is an integer. This relation can also be under- 
stood simply in terms of the round-trip condition in a 
cavity. For constructive interference to occur, the phase 
accumulated by the electric field in a round-trip [in this 
case 2{2Luon/c) -\- tt^ tt being the extra phase due to re- 
fiection from the two external mirrors] must be an integer 
multiple of 27r. When \Et^\'^ is plotted as a function of 
cj, it will then exhi bit p eaks with magnitude one at u 
values given by Eq. (A7). The line- width of these peaks 
is determined by the refiectivity of the external mirrors: 
the higher the refiectivity the narrower the peaks. 

What we have described up to this point is just a stan- 
dard Fabry-Perot cavity. If we now increase the refiec- 
tivity of the central mirror slightly, i.e. for |ro| <C 1, an 
electric field propagating along the x direction will also 
build up. To understand the physics of this system, one 
can think of the x and y propagating modes being weakly 
coupled by the central mirror. Then the eigenmodes will 
be symmetric and antisymmetric superpositions of these 
two modes, one along x and one along y as depicted in 
Fig. [6|a). The frequency separation between symmet- 
ric and antisymmetric combinations is proportional to 
the amplitude of the small enough refiectivity |ro|, which 
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can now be thought as the couphng coefficient, 
(a) 





FIG. 6: (Color online) Solid and dashed double- headed arrows 
symbolize two different modes that are weakly coupled by the 
central mirror, (a) |ro| ^ 1, (b) to ^ 1. 



The preceding ideas were exemplified by a numeri- 
cal calculation for which the total transmitted intensity 
\Et^ P + \Et^\^ + \Et^\^ as a function of uo is plotted in Fig. 
[Tlwith L — iTc/^^nuj^^ where uj^ is a convenient reference 
frequency. If Vo were identically zero, transmission would 
show peaks at frequencies which are odd integer multi- 
ples of cjr [c.f. Eq. (A7)]. For a small reflectivity, how- 
ever, one can notice that there occur two peaks around 
each of these frequencies. To verify that these two peaks 
correspond to the symmetric and antisymmetric super- 
positions of the two modes of Fig. [oj^a), we inspected 
the region around uo/uo^ = 1 more closely and plotted 
the magnitude and phase of the electric fields inside the 
cavity in Fig. [8) 
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FIG. 7: (Color online) Total transmission \Et\^ = \Et^\^ 4 
a function of uj/ujr with t = 0.1 and L - 



TTc/AnUr. Red dashed line shows the case for |ro| 
blue solid line is for Irol = 0.05. 



and 



One observation is that the peak positions and mag- 
nitudes at these peaks are nearly the same for all elec- 
tric field amplitudes. Nevertheless, one can identify two 
groups of fields in respect of their behavior as uj varies. 
The first group consists of £^i, £^5, and Eq contribut- 
ing to the y mode of Fig. [6) and the second group is 
composed of E'3, E4, £7, and E^ building up the x mode. 
One can notice from Fig. [S] that while the phases of a 
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FIG. 8: (Color online) Magnitude and phase of the complex 
valued electric fields inside the cavity as a function of cj/cjr 
for t — 0.1, |ro| = 0.05, and L — nc/AnuJr. Phases are defined 
in the range [— 7r,7r). y axis can also be used to read off \E\. 
Vertical dashed lines correspond to peak positions in the total 
transmission of Fig. ItI (a) Arguments of Ei^ E2, £5, and Eq 
are shown. Only |£^ijis depicted as other magnitudes behave 
very similarly, (b) Arguments of E3, E4, E7, and Eg along 
with I ^3 1 are shown. 



field in the first group are nearly the same at both peaks, 
the phase of a field in the second group at the higher 
frequency peak is shifted almost by tt with respect to its 
value at the lower frequency peak. This shows that the 
higher (lower) frequency peak corresponds to an antisym- 
metric (symmetric) combination of the x, y modes. 

Similar arguments can be advanced for the case to <C 1, 
i.e. when the central mirror is weakly transmitting. We 
can again describe the system in a two-mode approxima- 
tion and regard to as the coupling coefficient. As shown 
in Fig. |6|b), in this case the modes are confined in the 
lower and upper triangular regions of the cavity sepa- 
rated by the central mirror. This time peaks will appear 
around uj = 7rNc/2nL as opposed to the value given by 



Eq. (A7). This change is due to the phase acquired upon 
reflection from the central mirror. 

Two final remarks are in order, 
(i) As our scheme depends on the preservation of the 
type of circular polarization (defined in terms of angular 
momentum) inside the cavities, the central mirror should 
be chosen in a way that it reflects the two linear TE and 
TM polarizations almost equally and without introduc- 
ing any appreciable relative phase: this condition may 
be non-trivial to obtain for light at an oblique incidence. 
Still, any residual mixing of the (j± polarized eigenmodes 
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of the cavity can be suppressed by lifting their degener- 
acy, e.g. by introducing a slab of circularly birefringent 
medium in the cavity. 

(ii) Although we have restricted ourselves to a single cav- 
ity case in this Appendix, the same formalism can be 



generalized to several coupled cavities. As shown in [16], 
this leads to an alternative way of characterizing the op- 
tical response of a many-cavity system which does not 
rely on a tight-binding Hamiltonian such as 
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